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Abstract 

We investigate the contributions of quantum fields to black hole 
entropy by using a cutoff scale at which the theory is described with 
a Wilsonian effective action. For both free and interacting fields, the 
total black hole entropy can be partitioned into a contribution derived 
from the gravitational effective action and a contribution from quan- 
tum fluctuations below the cutoff scale. In general the latter includes 
a quantum contribution to the Noether charge. We analyze whether 
it is appropriate to identify the rest with horizon entanglement en- 
tropy, and find several complications for this interpretation, which are 
especially problematic for interacting fields. 

1 Introduction 

The concept of black hole entropy, first introduced by Bekenstein [H |2] us- 
ing information theory and the analogy between black hole mechanics and 
thermodynamics, became firmly established when Hawking [3] derived black 
hole radiation and its thermal properties by considering quantum fields on 
a black hole background. The formula expressing the black hole entropy in 
terms of the horizon area A and Newton's constant. 
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has since then become a focal point of quantum gravity research. Attempts 
to derive it from an analysis of microscopic, fundamental degrees of freedom 
have met with varying degrees of success in string theory, loop quantum 
gravity, and other approaches. 

Staying within an effective treatment in which gravity is described by 
a metric field, there appear to be two contributions to black hole entropy. 
Firstly, the gravitational field itself in the absence of matter fields seems 
to have a "gravitational" entropy given by the Bekenstein-Hawking formula. 
This can be derived from the first law of black hole mechanics and the Hawk- 
ing temperature, or from the saddle-point evaluation of the Euclidean grav- 
itational path integral for the thermal partition functionjl]. Secondly, the 
contribution of quantum matter fields on the black hole background (and of 
gravitons, if metric fluctuations are quantized perturbatively) to the thermal 
entropy is also proportional, in the leading order, to the event horizon area, 
but with a divergent coefficient. This entropy arises from the one- loop cor- 
rection to the thermal partition function. It can be conceptualized as the 
entropy of a thermal state for quantum fields outside the horizon, or, at least 
in some cases, as the entanglement entropy across the horizon of quantum 
fields in a global pure state. For reviews of the main issues and results, see 

The area-scaled divergence in the matter contribution to the entropy can 
be absorbed into a renormalization of Newton's constant. More precisely, the 
divergences in the entropy are related to the "bare gravitational" entropy in 
the same way as the divergences in the effective gravitational action are 
related to the bare gravitational action [7]. The renormalization properties 
of the black hole entanglement entropy have been further studied in [HI [9|[T0]. 
with the case of non-minimally coupled fields receiving special attention both 
for scalars [TT] and for gauge fields [121 [131 dl]- Different regularization 
methods (brick-wall boundary at the horizon, Pauli-Villars regulator, UV 
cutoff e in heat kernel expansion...) do not always give the same results for 
the matter contribution [H [151 [HI 113 EHl ESj- These discrepancies are not 
surprising, since the UV regulator modifies precisely those degrees of freedom 
which are most responsible for the entropy. 

While black hole entropy emerges, formally, from the gravitational parti- 
tion function, that does not reveal the nature of the states that are counted 
by the entropy. It is tempting to think that, like the minimally coupled 
matter contribution, all of the entropy might be interpreted as entanglement 
entropy. For this to make sense, it would seem that the low energy Newton 
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constant must arise fully from integrating out quantum fluctuations, i.e. that 
there is no "bare" gravitational action or entropy at the UV cutoff scale, so 
that gravity is entirely "induced" [10] . From the QFT viewpoint there seems 
to be no reason why this should be so, although a thermodynamic argument 
suggests it must be [20], [21]. But, in any case, without a UV completion of 
the theory, it is not really possible to assess this entanglement interpretation 
of the full entropy because the value of the entropy depends on the artificial 
UV cutoff at scale A^/y. 

However, it is possible to test the entanglement interpretation of black 
hole entropy in a way that side-steps the unkown physics of the UV cutoff and 
deals only with finite quantities. We can partition the degrees of freedom into 
those with momenta greater than an intermediate scale k <C A[/y, which are 
integrated out and absorbed into a Wilsonian effective action, and those with 
momenta less than k, whose quantum fluctuations contribute explicitly to the 
entropy. As k is lowered, contributions to the total entropy transfer from the 
explicit fluctuations to the "gravitational entropy" of the effective action, i.e. 
the area term plus curvature corrections, via the flow of the gravitational 
couplings. This allows us to exhibit the renormalization behavior of the 
contributions to the entropy without sensitivity to the UV cutoff, which in 
turn allows us to make precise sense of the question whether the entropy of the 
modes with momenta less than k admits an interpretation as entanglement 
entropy. We shall find that such an interpretation, while at first superficially 
plausible, suffers from a number of difficulties. 

The renormalization group (RG) flow of black hole entropy from state 
counting of explicit field fluctuations to effective gravitational action was 
studied long ago in Ref. [22]. That paper studied how the entropy account- 
ing changes when the renormalization scale drops below a mass scale of the 
fluctuations, and focused on the interpretation of the contribution from non- 
minimal couplings to curvature. Our study is very similar in spirit, but we 
consider continuously varying RG scale, and we try to assess more precisely 
the validity of the state counting interpretation even when non-minimal cou- 
pling is absent. 

A somewhat similar scheme was introduced recently in [23] (see especially 
subsections 4.1.3 and 4.3.3). The main difference is that in that paper the 
full physical description of the system is assumed to be contained in the 
gravitational effective action, with no consideration of the quantum fluctu- 
ations below the cutoff scale. The same assumption is made in [S^, where 
an explicit identification of the cutoff scale with the size of the black hole 
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is proposed. We take an alternative interpretation of the renormalization 
group, in which the cutoff scale is an arbitrary parameter separating short 
and long wavelength modes, and use it as a tool for probing the entanglement 
interpretation of the long wavelength contribution to the entropy. 

The structure of the paper is as follows: In Section 2 we will present the 
definition of black hole entropy from a canonical ensemble and its connection 
with the effective action for the gravitational field. In Section 3 we focus on 
free quantum fields and we introduce the definition of the running effective 
action dependent on a RG cutoff scale, and exhibit the renormalization prop- 
erties of black hole entropy, including a detailed computation for the massless 
scalar field. In Section 4 we present a scheme for extending this idea to in- 
teracting fields. In Section 5 we discuss a number of issues that impede a 
direct interpretation of the IR contributions to the entropy as entanglement 
entropy. Section 6 includes a summary of the results and a discussion. 

Throughout this paper we work in four dimensions and use units with 
h = c=l. 

2 Canonical ensemble and entanglement en- 
tropy 

This paper will focus on the properties of black hole entropy in a thermal 
state, as defined by a canonical partition function at fixed temperature: 



The canonical Hamiltonian "H in (|2]) includes terms both for the gravita- 
tional field and the matter field. Because of diffeomorphism invariance, l-L 
is a boundary term, when acting on physical states that satisfy the diffeo- 
morphism constraints. The thermal ensemble is thus specified by boundary 
conditions. The entropy is computed, as is standard in statistical mechanics, 
by application of the operator {Pdp — 1) to — InZ. This expression for Z is 
formally equivalent [1] to the Euclidean path integral 



Z = Tr e-^^ . 



(2) 




(3) 



4 



where Sf, is a bare action for the gravitational field and S{g, (p) is an action 
for the matter field if on background The notation /^^^ represents in- 
tegration over Euclidean fields with the metric gs fixed and stationary with 
periodicity /3 in Euclidean time at an outer boundary. For the matter fields, 
we are free to choose any boundary condition (e.g. Dirichlet or Neumann) as 
long as it is stationary and compatible with the /3-periodicity. These alter- 
native choices may represent genuinely different physical ensembles leading 
to different results for the entropy. In this paper we will make the simplest 
assumption of Dirichlet boundary conditions, leaving a broader discussion of 
the issue for later work. We take the boundary to have a finite size, small 
enough for the canonical ensemble to be stable [25]. (Alternatively we could 
work in asymptotically Anti-de Sitter spacetime[27].) 

We can formally integrate out the matter field in ([3]), defining 

W[g] = - In / Vcp e-^t^''^] , T[g] = St,[g] + W[g] , (4) 

so that 

Z{f3,gB)= [ Vge^^^^^. (5) 

The matter contribution W to the gravitational effective action F is generally 
UV divergent. We assume that a regularization scheme for it is in place, and 
that the bare gravitational couplings (parameters of Sb) are adjusted so that 
the renormalized couplings in F are finite. 

We use now a zeroth-order approximation for quantum gravity, in which 
the gravitational path integral is evaluated at the saddle point. (This 
amounts to disregarding from the calculation all graviton fiuctuations; in 
a more full treatment, which we omit for simplicity, they could be included 
perturbatively among the matter fields ip.) Therefore we write: 

^(/3,^?B)=e-r[5(^'^-)]. (6) 

Here g{P,gB) is the metric that solves the equations of motion derived from 
T[g], with the boundary conditions that we have set. We assume that the 
renormalized cosmological constant is zero or negligible, and that the higher- 
order in curvature terms of F (which include nonlocal terms) can likewise be 

■^A demonstration of this formal equivalence taking the constraints into account is given 
in Ref. [21]. 
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neglected in the regime of interest. Then r[g] is comprised only of the bulk 
curvature term and the corresponding Gibbons- Hawking boundary term0: 

Now we specialize the discussion to the ensemble defined by a spherically 
symmetric boundary metric which is a 2-sphere of radius r^. We also 
work in four dimensions, though most of our results generalize straightfor- 
wardly to d dimensions. The on-shell metric g is Euclidean Schwarzschild: 



- \ -1 



ds^=(l_!±)dt2+(l_!±) dr^ + rMl]^, (8) 



with the horizon radius f+(/3,rB) defined so that (i) there is no conical sin- 
gularity at r = r+ and (ii) the ensemble is stable ((ir+/(i/3 < 0). These 
conditions imply [SB] that r+ is the larger root of the equation: 

/3 = 47rr+(^l-^y . (9) 

The "box" at the boundary tb stabilizes the ensemble by giving the black 
hole a positive heat capacitylfl 

When evaluated on this metric, the bulk term vanishes and the boundary 
term yields: 

r[^] = T^(37rr^-47rrBr+), (10) 

We omit the latter because it is proportional to /3 and hence does not affect 
the entropy, which is given by the standard thermodynamical formula 

Sbh = - {(3dp - 1) In tb) = {f3dp - 1) r[^(/3, r^)] . (11) 



^It has recently been pointed out ^23j that the Gicn that appears in the boundary term 
of the effective action might not be the same as the one in the bulk term when the matter 
field is non-minimally coupled to the curvature. For the moment, we assume that this is 
not the case and that bulk and boundary Gron are equal, as is the case, for example, when 
Lp is the minimally coupled scalar field. 

^If the black hole grows a little bit, its horizon is closer to the box, so less redshifting 
of temperature occurs from the horizon to the box. If < 3r+/2, then this suppression 
of redshifting makes the temperature of the larger black hole higher than 1 / (3 at the box. 
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When applied to fllUI) . this results in the renormalized Bekenstein-Hawking 
formula: 

Sbh = j^. (12) 

where A = 47rf^ is the horizon area. 

The above approach to evaluating the matter field contribution to Sbh 
has been dubbed in the past the "on-shell" computation or the "thermo- 
dynamical" computation [161 HI]- Within this procedure, the fact that the 
full entropy including quantum corrections is expressed by the Bekenstein- 
Hawking formula involving the renormalized Newton constant Gren is an 
immediate consequence of the renormalization of the effective action, as was 
emphasized in [H]. Note that this also implies that the so-called "species prob- 
lem" (the dependence of the quantum contribution to black hole entropy on 
the kinds of existing quantum fields, apparently contradicting the universal- 
ity of the Bekenstein-Hawking formula [28]) is moot because in terms of the 
renormalized value Gren, this formula is always correct regardless of the field 
species (which affect the relation of Grcn to the unobservable bare value). 

In contrast, when computing the contribution of the matter fields as en- 
tanglement entropy on the black hole background [281 ISll |30] , this renormal- 
ization property is much less apparent. We shall therefore discuss in more 
detail the relation to entanglement entropy computations. 

To compute the entanglement entropy contribution to black hole entropy, 
we consider a minimally coupled quantum field on a Schwarzschild back- 
ground. The entanglement entropy across the event horizon is defined by 

Sent = -Tr Pout In Pout , (13) 

where pout is the restricted density matrix for the external region, with the 
internal states traced over. Using a Euclidean path integral representation 
for p, this can be rewritten [211 E21 123] as the operation {ada — 1) applied 
to the matter contribution Wa to the effective action, computed on the Eu- 
clidean Schwarzschild background with a conical singularity introduced at 
the horizon. Here 27ra is the periodicity of the angular "time" coordinate 
and a is set to 1 after the differentiation: 

(14) 

a=l 

Actually this "conical" procedure yields the entanglement entropy only for 
minimally coupled fields. For non-minimally coupled fields, Scone is not equal 



5'cone = («5q - 1) 
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to S'cnt, as it includes an extra, contribution interpretable [Ml [13] as the ex- 
pectation value of a term in the Wald entropy [53] . Our subsequent discussion 
in this section assumes minimal coupling. We shall return to expand a bit 
on this issue in Section 5. 

The use of a background with a deficit angle 27r(l — a) is equivalent to 
using periodic quantum fields whose physical periodicity /3 at a radius is 
related to the horizon parameter r_,_ by 

/ \ 1/2 

/5 = 47rr+a f 1 - , (15) 

rather than by (Q. Hence a and /3 are proportional at fixed r_|_, and the 
entanglement entropy is equivalently given by 

(16) 

r+=r+ 

Comparing with (ITT!) and recalling the relationship between the en- 
tanglement black hole entropy and the thermodynamical black hole entropy 
becomes clear: The entanglement entropy (when computed using as back- 
ground the metric g) differs from the matter contribution to the thermody- 
namical entropy, only in that the variation with respect to /3 is done keeping 
r_|_ fixed (thus involving the introduction of a deficit angle), instead of taking 
into account the dependence of r+ on /3 in the on-shell solution Hence the 
dubbing of this as the "off-shell" method in [T6] . 

However, it can be argued that the two methods give actually the same 
results. First, consider again f lTT]) . and recast the dependence of T[g{(3,rB)] 
as T[g{f3,rB,f+{f3,rB))]- Then the "total" /3-derivative appearing in (fTTj) can 
be unpacked as 

dpT[g{/3,rB)] = (^d^T[g{f3,rB,r+)] + dr^T[g{(3,rB,r+)]dpr+^ 

The last term should vanish, because dr^T[g{P,rB,r^] expresses a variation 
of the action with respect to the metric with fixed boundary conditions, 
which is zero at the on-shell value of r_|_ (a stationary point of vj^. Indeed, 
this can be checked with a direct computation of this derivative evaluated 

^This argument has been made previously in [llj : see also [36] . 



S,^, = {f3dt,-l)W[g{/3,rB,r+)] 



(17) 
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on the deficit-angle version of Scliwarzscliild, wliere tlie conical singularity is 
accounted for as an extra term in the curvature scalar: 



R = R + 4:7i{l-a)6j:. 



(18) 



So we conclude that the thermodynamical entropy can be expressed using 
just the first term of (fT71) . which includes the same kind of "off-shell" partial 
derivative as (ITB]) . Hence we have: 



Therefore the entanglement entropy must share the renormalization property 
of the thermodynamical entropy: its divergences should be absorbable in a 
redefinition of the bare gravitational couplings, as it occurs for W in the 
thermodynamical entropy. 

Note that the equality requires the entanglement entropy to be com- 
puted using as background spacetime the metric g that solves the quantum- 
corrected equations of motion. In practice, this metric is assumed to be 
Schwarzschild (or another known black hole solution to the gravitational 
theory) expressed in terms of the observable, low-energy couplings. 

An issue we have mentioned but not dedicated proper attention to is the 
need to regularize W (the matter contribution to F) to make it finite. Insofar 
as W is divergent, the whole argument is not rigorously defined. It would 
clearly be preferable if the renormalization properties of the entropy and 
the relation to entanglement entropy could be studied by manipulation of 
manifestly finite quantities only. In the rest of this paper, we achieve this by 
introducing a Wilsonian renormalization group scale k, at which the entropy 
can be decomposed into respective contributions from the gravitational action 
and the matter action. The flow of these contributions as k changes is then 
well-defined. In the next section we study this flow for free matter fields, 
and in section 4 matter interactions are included. In section 5 the viability 
of interpreting the matter contributions as entanglement entropy is probed. 

3 RG flow of black hole entropy for free flelds 

In this section we will show how the black hole entropy, in the "thermody- 
namical" framework presented above, can be described in a way that makes 



Sbh = {Pdp - 1) {S,W, rs, r_,)] + WW, rj,, r+)]) 
= (/39^-l)5,[(7(/3,rB,r+)] _ + S-ent • 



(19) 



r+=r+ 
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its renormalization properties clear avoiding the handhng of divergences. In 
the spirit of the Wilsonian interpretation of the renormahzation group [H7] . 
the idea is to introduce a cutoff scale k and to integrate out only the quantum 
modes above that scale. 

Let us start by going back to (Q, after the zeroth-order approximation 
for the gravitational path integral has been made, and replace r[g] by its 
definition: 

= e-5,[3(/3,rs)] f J)^e-S[mrB)M _ (20) 
Jl^,rB 

The path integral should be assumed to contain an implicit covariant UV 
regulator, for example a short-distance cutoff e in the heat kernel expansion. 
If <y9 is a free field, assumed for illustration to be a scalar, then its action 
takes the general form0 

S[g,^] = J y/g^^g^, (21) 

where is an elliptic operator depending on the metric g; for example A = 
— for a massless minimally coupled scalar field, and A = —V^+^R{g) for a 
nonminimally coupled one. Here is the Laplacian operator on background 
g. The matter contribution W to the effective action is given by the one-loop 
determinant: 

T[g] = S,[g] + ^TrJnAg, (22) 

which can be computed with standard heat kernel expansion techniques. 
With the regulator e removed, the right hand side is divergent; leaving it in 
place, we can make the couplings of the bare action St dependent on e in 
such a way that the effective action F and the physics derived from it are 
independent of e. 

We now want to introduce an intermediate RG scale k < 1/e. One simple 
way of doing this is using an additive cutoff as introduced by Wetterich [3H] 
to study the exact renormalization group. The gravitational effective action 
at scale k is defined in general as: 




^We assume that there is no boundary term in the matter action, or more exactly that 
it vanishes when the boundary conditions are imposed. 
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Here7^fc(A) = k^r{§), with r{z) being a function that satisfies the properties 
r(0) = 1 and r{z) =0 for z >1. This imphes that the TZk term serves as an 
IR cutoff in the path integral, suppressing from the contribution of the 
modes with eigenvalue p"^ < k"^. Hence is an "average" effective action that 
only incorporates the effect on the gravitational couplings of fluctuations on 
length scales smaller than k~^. For — j- 0, F^ approaches the full effective 
action 

In the free fleld case, F^ can be computed exactly as a modifled one-loop 
determinant: 

^k[9] = St[g] + ^Tr, In [A^ + 7^fc (A,)] , (24) 

which by differentiation with respect to RG flow "time" t = Ink yields the 
well-known RG flow equation 



dtTk = -Tr 



dtTZ, 



Ay + 7^fc (Ag 



(25) 



Comparing flMl) and fl22l) . we get the following expression for the difference 
between the total effective action and the effective action at scale k: 



r[g]-T,[g] = ^TT\n 



A. 



Wk[g]. (26) 



Notice that we can drop the UV cutoff e in this expression, because the trace 
has now acquired a lower intrinsic UV cutoff at scale k since TZk = for 
eigenvalues > k"^; the right hand side of (126|) is well-deflned regardless of 
any overall UV cutoff for the theory. Evaluating at the on-shell Euclidean 
Schwarzschild metric g = g{P,rB) and applying the operator {(38/^ — 1) we 
obtain the relation for the entropies: 

Sbh-SPj,= (^f3dp-l)Wk. (27) 

Sbh is the total black hole entropy, incorporating the bare gravitational 

(k) 

contribution and the total effect of the quantum matter flelds. S"^^ is the 
"effective gravitational entropy at scale k" , derived from the effective action 
Ffc[5f]; it incorporates both the bare gravitational contribution and the ef- 
fect that the high-frequency quantum modes have on the renormalization of 



^Notc that with the cutoff function TZk the suppression of the modes with < is 
not complete; it is rather hke giving them a mass ^ . 
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the gravitational couplings. The right hand side encompasses the effect of 
the low-frequency modes, whose contributions had been suppressed from F^. 
The total entropy Sbh is independent of the sliding scale k, which partitions 
it into two complementary contributions. Thus the renormalization proper- 
ties of the entropy are made clear without need to worry about the global 
UV cutoff and the treatment of divergences. Sbh is expressed in terms of 
low-energy, observable couplings, and the other terms in the equation dif- 
fer from it by a finite RG scale k, which in principle can receive a physical 
interpretation. 

To investigate whether the right hand side of ( 127|) can be interpreted as 
the entanglement entropy of the modes below scale k, let us evaluate it in a 
concrete example. We consider a minimally coupled massless scalar field (p, so 
Ag = — V^. We choose as cutoff function the "optimized cutoff" introduced 
by Litim [39], which is given by: 




We have in this case that the operator traced upon in Wk (126|1 is exactly zero 
for eigenvalues above fc^: 



We can use the procedure described in the Appendix A of [IQ], adapted 
for a 4-dimensional manifold with boundaries, for using the heat kernel ex- 
pansion to compute the trace of an arbitrary function of our operator, F{A). 
Assuming Dirichlet or Neumann boundary conditions are imposed, one has 





(28) 




(29) 



Tr[F(A)] 



1 r 



ao{A)Q2{F) + ai(A)g3 (F) + a2{A)Q,{F) + ■■■ 



167r2 L 



(30) 



where the coefficients Qn{F) are defined as 




(31) 
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and the heat kernel coefficients take the form [IT] : 



ao(A) = 
ai(A) = 
a2(A) = 




For F{z) = In [p-] B [fc^ — 2], we have Qn = —k'^'^/'n?V{n). Hence we obtain: 



(35) 

The next step is to evaluate at the on-shell Euclidean Schwarzschild met- 
ric g{P,rB) and to compute the right hand side of (1271) . The bulk curvature 
term is zero, and the bulk and boundary volume terms are proportional to 
P and hence vanish upon application of the {Pdp — 1) operator. The bound- 
ary K term gives a result proportional to the event horizon area. Hence, 
neglecting higher-order curvature terms, f l27|) takes the form 

A A _A I' e 

This equation can be read in two ways. On one hand, canceling the y4/4's, it 
just expresses the RG running of G due to the quantum corrections induced 
by the scalar field, which was already implicit in fl35|) . On the other hand, 
we can interpret it as expressing two contributions to black hole entropy: 
For any scale k, the total black hole entropy A/4Go (where Gq is the fully 
renormalized Newton constant) can be partitioned in two contributions, the 
effective gravitational entropy at scale k (which is A/4Gk) and the contribu- 
tion of the scalar field's modes that are below k, which is A k'^/A87T. When 
sliding the RG scale k, the balance of the entropy is shifted between the two 
terms, leaving the total entropy unchanged. 

It would seem natural to regard the contribution of the lower modes as 
corresponding to their horizon entanglement entropy. It has the same form as 
the the total entanglement entropy calculated with a UV momentum cutoff 
A, with the intermediate scale k playing the role of A. Of course the precise 
expression for the running of Gk depends on the cutoff function TZ^. A 
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327r2 



-R 



dM 



Ak-" 
~9~ 



-K 




different choice of regulator would lead to a different numerical coefficient of 
Ak"^. In itself this does not seem problematic for the entanglement entropy 
interpretation, since it just reflects the implementation of the partitioning of 
the contributions from degrees of freedom above and below the RG scale. 

In section 5 we shall discuss some other questions concerning the justi- 
fication of the entanglement interpretation of the contribution of the lower 
modes. First, however, we consider how the preceding analysis must be mod- 
ified in order to account for interactions of the matter degrees of freedom. 

4 Interacting fields 

In this section we will discuss how our framework can be extended to inter- 
acting quantum fields. It turns out that the distinction at a scale k between 
the gravitational entropy and the contribution from the modes below k can 
be defined as we did for free fields, though it is much more difficult to write 
down the exact form of each contribution for a given example. There are 
also further interpretational issues, which will be addressed in Section 5. 

Let S[g,(f] be a bare action for the quantum field (p on the gravitational 
background g. Equation (I2UI) for the black hole partition function is true 
regardless of whether S contains interactions. We can, as before, introduce 
an additive cutoff TZk and define the running gravitational effective action 
^k[g] by f l23|) . and define the gravitational part of the entropy at scale k by 
applying the (/39/3 — 1) operator to rfc[^]. However, since for interacting ip the 
one-loop determinant is not an accurate evaluation of the effective action, we 
are lacking a compact expression like fl27l) for the contribution of the modes 
below k. Integrating out the upper modes not only produces running in the 
gravitational effective action, but also the Wilsonian effective action for ip 
depends on the scale k. In the following we elaborate on how to quantify 
this running and obtain expressions for interacting fields as close as possible 
in spirit to (1271) . 

Let us start, again, with the total partition function, with an overall UV 
cutoff implicitly in place with a short- distance regulator e. We isolate the 
kinetic term in the matter action and define as Sb[g, ^p] the sum of the bare 
gravitational action Sii[g] and the non-kinetic terms of the bare matter action. 
Then the partition function is given by 

Z[g] = e-^[^?] = /l^cpe-^^^^-^^V-^i^b^v] _ (37) 
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We now introduce the IR cutoff function 71^, witli the same properties as in 
the previous section, defining the Wilsonian effective action at scale k by: 



-Sklg,<t>] 



Dip e 



y ip{-\'l+Tik{-\'l))v-Sb[g,'l>+v] 



(38) 



The purely gravitational part of Sk coincides with our previous definition of 
Ffc, and we shall be decomposing the entropy into the contribution from 
and that from the remainder, which captures the physics of the lower modes 
of the matter field (p on the background g. To that end, we introduce the 
notation 



^ Tklg] 
Sk[g,(p] 



Sk[g, 
Sk[g, 



0] 

-Tklg]- 



(39) 
(40) 



The action Sk includes the effects of the modes "above k" , so the partition 
function (1371) should be expressible as a path integral over the modes "below 
k" using this action. We find such an expression by assuming it can be 
written in the form 



Z[g] = N,[g] / P0e-^ 



.i/0p-V-5fc[9,' 



(41) 



for some suitable choice of normalization Nk and IR propagator P^. Substi- 
tuting f p8|) into f HTj) . shifting one of the field variables so that its integral is 
Gaussian, and performing the integral, we find that fl37|l is recovered with 
the following definitions of and P^: 



Nk[g] = det^/2 



(-v^ + 7^,(vp)^ 



^ 9 



-v^ + 7^,(v2) 



(42) 



(43) 



Note that since TZk vanishes for modes above k, so does the propagator P^ 
and hence the path integral in (HTl) acquires a UV cutoff at this scale. We 
stress that ( HTjl is identical to the full path integral, with the information 
about the upper modes encoded in the Wilsonian effective action Sk[g,(l)], 
which is obtained from integrating them out according to fl38|) . 

The definition of Sk given by ( l38l) is purely formal however, and unsuit- 
able for analyzing the entropy of the black hole at scale k. In first place, the 
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path integral cannot be computed in a closed form for an interacting the- 
ory. Moreover, the expression in terms of a bare action and a path integral 
which is divergent requires an exphcit regularization procedure to deal with 
divergences. This goes against the spirit of our approach, based on analyzing 
the difference between the expressions for the entropy obtained at different 
effective scales, in terms of finite, effective quantities only. 

The right tool for these purposes is an RG flow equation, detailing how 
Sk[g, 0] changes with the scale /c in a local way. This is the Polchinski equa- 
tion [42], which in the present case takes the form: 




SSk 



Tr 



Pk- 



Tr 




where the overdots represent ^-derivatives, and the center dot notation (■) is 
explained in Appendix A. Since this equation differs by the last term from 
the fiat space form that is derived in standard presentations of the Wilsonian 
renormalization group |l3l HI] , we detail in Appendix A how it is obtained 
using (138|) as starting point. 

Expanding the flow equation with a systematic approximation method 
(e.g. a derivative expansion) would give beta functions for each of the grav- 
itational couplings in Tk and of the field couplings in Sk- Let us assume 
we are in possession of a solution to these flow equations for the couplings 
(found, perhaps, with numerical techniques, and using as initial condition for 
the flow a known form of the effective action at low energies). This would 
then allow us to write down the form of Tk[g\ and Sk[gi 0] for any given value 
of the RG scale k. The log of the partition function ( 14T]) is then expressible 
in terms of these quantities as 



\r.Z[g] = T[g]=Tk[g] + Wk[gl 



with 



Wk[g] 



In 



(46) 
(47) 



In this way the free energy is decomposed into finite, purely gravitational and 
matter parts in a scale-dependent manner. This decomposition generalizes 
the one (126|1 discussed previously for free fields. 



16 



The saddle point approximation to the total entropy of the thermal 
ensemble is obtained using the partition function evaluated at the metric 
g{P, tb) given by ([8]), which is a solution to the full effective action at k = 0, 
i.e. with all fluctuations integrated out. Using the decomposition P6 | Wf} for 
Z, we obtain: 

SBH = {l3dp-img] + Wk[g]) 

=: Sl,% + iPd^ - l)W,[-g] ■ (48) 

The term S^^jj, is (as before) what we define as the gravitational black hole 
entropy at scale k. The second term, with the definition ( 147|1 . generalizes 
the right hand side of ( 1271) with its definition (126|) . and encompasses the 
contribution to the entropy of the quantum modes below scale k, as computed 
with the appropriate Wilsonian action for them. It can easily be checked that 
for the case of free fields, where the path integral is Gaussian and can be done 
exactly, f H7|) reduces precisely to f l26|) so that both expressions for the entropy 
agree. 

Once more, the renormalization properties of the entropy are made ex- 
plicit without the need of specifying the global UV cutoff and handling di- 
vergences; a global implicit regulator e is needed to make flTTI) well-defined, 
but it drops from the calculations and when we reach (148|) we are dealing 
only with finite physical quantities. The total entropy, the right hand side 
of (l48l) . does not depend on the scale k, but when the scale is shifted the 
relative balance of the two contributions on the right hand side is changed, 
as the gravitational and non-gravitational part of the effective action flow 
according to (HSj) . 

5 Interpretation of the IR contribution to the 
entropy 

We have studied in the previous two sections the RG flow of black hole 
entropy contributions coming from above and below a running scale k. We 
obtained expressions fl27|) for free fields and fHHl) for interacting fields, with 
the associated definitions fl26llTn) for Wk in each case. These decompositions 
are interesting in their own right, as they illustrate how a black hole entropy 
computation tracks the RG flow of gravitational couplings. But our main 
motivation for undertaking this exercise was to test the notion that black 
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hole entropy is, at least in part, entanglement entropy of quantum fields, in 
a controlled setting where no divergent quantities arise and where properties 
of the UV completion of the theory are irrelevant. In this framework, the 
results may have a more direct and less ambiguous physical interpretation. 

In particular, the tempting interpretation of fl27j) and pHj) (as discussed 
briefly with regard to fl5B|) ). is that the IR contribution to the entropy. 



can be identified with the entanglement entropy of the modes below scale 
k. There are a number of considerations that complicate this interpretation 
however. We will now discuss them. 

5.1 Contact terms 

As mentioned in section 2, for non-minimally coupled scalar fields, and per- 
haps for gauge fields and gravitons, the "conical entropy" (i.e. the entropy 
of the thermal partition function defined on a space that acquires a conical 
deficit when the Euclidean period is varied off-shell) contains a contribution 
from the tip of the cone, the so-called contact term, that does not appear 
to admit a statistical interpretation]^ This contact term arises also in the 
contribution from the modes below scale k, so in general that contribution 
would not consist only of entanglement entropy. In a specific interacting 
model in 1-1-1 dimensions [22], it was illustrated how non- minimal coupling 
and the associated contact term can arise from a Wilsonian effective action 
when some degrees of freedom above a certain mass scale are integrated out. 
Hence it is possible, in a given setting, that the contact term is a stand-in 
for an entanglement contribution, but that need not be so. 

The contact term is a hybrid between a gravitational and a quantum 
contribution. It can be interpreted [SU [13] as a term in the Noether charge, 
i.e. in the Wald entropy |3S], involving the expectation value of the squared 
matter field with the cutoff k. Thus, in the presence of non-minimal cou- 
pling, we should refine our conjecture about the decomposition of the total 
entropy at scale k into gravitational and entanglement contributions. The 
gravitational part must include all contributions from the Noether potential 
at scale k. For example, for the scalar field with an i?</?^ coupling in Sk, the 

^However, in Refs. [3411111 US] a statistical interpretation was proposed in terms of zero 
energy modes localized at the horizon. 




(49) 
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gravitational part would arise both from this term and from r^. Since (f is 
a fluctuating quantum field, this makes the Noether potential an operator 
rather than a classical quantity. The gravitational part in the conjecture 
involves the expectation value of this operator. 

5.2 Euclidean vs. Lorentzian RG scale 

The running scale k in our calculations is defined as a cutoff in the eigen- 
values of the Euclidean Laplace operator. In the Euclidean domain, we have 
a clean separation between two contributions to the entropy: the term of 
(HSj) coming from Ffc represents the "gravitational entropy" at the scale k, 
and the one involving 5*^ represents the contribution of the Euclidean modes 
below k. But what exactly do these terms correspond to in the Lorentzian 
domain, where the entanglement entropy is fundamentally defined? It seems 
there is a close relationship, but it has not been spelled out0. Until the 
relationship between the Euclidean and the Lorentzian domains is clarified, 
it is questionable to identify the quantity f l49l) . defined by a Euclidean cut- 
off, with the entanglement entropy of a subset of the Lorentzian quantum 
fluctuations. 

5.3 Uncertainty relation between horizon location and 
momentum cutoff 

The notion of "horizon entanglement entropy" refers to the von Neumann 
entropy of the reduced density matrix of the exterior degrees of freedom. 
We are here considering this notion in the presence of a momentum cutoff. 
If the calculation were strictly in the Lorentzian domain, the limitation on 
momenta would presumably imply that the separation of degrees of freedom 
on one side of the horizon could not be arbitrarily sharp, and in fact would be 
fuzzy at the scale of the inverse momentum cutoff. Since the entanglement 
entropy is dominated by the contributions at the shortest scales, this means 
that there would be an order unity fuzziness in its value, computed in this 
fashion. This dependence on the cutoff would not be so disturbing if it were a 
feature of a UV regulation of an otherwise undefined quantity, but one might 

"'^'^In connection with this, one interesting avenue for future research is to follow the 
methods of [47] , where the RG cutoff is defined in terms of the eigenvalues of the spatial 
part of the Lorentzian Laplace operator in a 3+1 decomposition. 
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have hoped for a more precise definition of the entanglement entropy at the 
RG scale k. 

We do not encounter this fuzziness in our computation. The unregulated 
Euclidean path integral for the partition function can be viewed as a formal 
computation of the trace of exp{—PHext), where He^t is the Hamiltonian 
for the degrees of freedom exterior to the horizon. However when this path 
integral is filtered to include only the modes with momentum below k, it 
is no longer exactly the trace of an operator on the exterior Hilbert space. 
For this reason, what we compute in (1491) is not precisely the entropy of a 
reduced density matrix corresponding to a subset of the exterior modes. 

5.4 Momentum entanglement 

For interacting theories, a further complication arises. In the ground state, 
degrees of freedom with different momenta are entangled. In Minkowski 
spacetime, the reduced vacuum density matrix for IR degrees of freedom be- 
low a scale k has an entanglement entropy per unit volume. This was com- 
puted perturbatively for various theories in various dimensions in Ref. |1H], 
but it is explained there that for some theories the perturbative calculation is 
not adequate. Nevertheless, we can estimate that a lower bound for theories 
in 3+1 dimensions should scale as A^/c^, where A is the coupling constant. 
The actual result for a given theory might involve some power of a higher 
energy scale M and logarithms of the ratio M/k. 

Momentum entanglement will also play a role for the reduced density 
matrix of IR modes outside a black hole horizon. If we restrict attention 
to the volume of space at a proper distance / from a spherical horizon, the 
lower bound for momentum entanglement entropy would scale as X^k^lr"^, 
whereas the horizon entanglement entropy scales as k^r"^, hence the former 
dominates unless X'^kl < 1. If A is much smaller than unity, the momentum 
entanglement contribution could potentially be suppressed by focusing only 
on a region of radial width / < (A^A;)^^, which would be much larger than 
the cutoff wavelength and hence compatible with the cutoff. 

The IR contribution (H^ to the total entropy does not appear to have 
any contribution corresponding to momentum entanglement, and it should 
not, since that is "internal" entanglement that does not contribute to the 
total entropy. Therefore, fH9|) must differ from the von Neumann entropy of 
the reduced density matrix of the lower modes outside the horizon. We now 
attempt to identify the origin of this discrepancy using a formal computation 



20 



in which the issues raised in the previous two sections are ignored. This strat- 
egy is sensible because the issue of momentum entanglement is orthogonal 
to the others. 

Let Z = Tra,A exp(— ) denote the full partition function of the exterior 
degrees of freedom, where a and A stand for IR and UV degrees of freedom, 
schematically. If we first trace only over A, we have 

TvAeM-PH) = Zgexp{-PHa) . (50) 

Here is a /3 dependent number, independent of the fields a, and Ha is an 
effective Hamiltonian for the lower modes. This split may be ambiguous in 
general, but for the scalar field we defined Zg = exp{—rk[g]) via the part of 
the effective action that was independent of the scalar field, and the remaining 
effective action was Sk, which is the action that would correspond to Ha- In 
the presence of interactions we expect Sk to include non-local terms, and we 
expect (3 dependence in Ha simply because it is defined by a /3 dependent 
procedure. 

It follows then that Z = ZgZa, where Za = Traexp{—/3Ha). Now when 
we compute the contribution to the entropy, —{I3dp — l){]xiZg + InZa), the 
Zg term contributes a "gravitational entropy" and the Za contribution corre- 
sponds to din])- If Ha did not depend on /3 (for example for a non-interacting 
field), then the Za term would contribute the von Neumann entropy of the 
density matrix pa = Z'^ exp{—f3Ha). This would just be the entropy of 
the a subsystem. However the /3 dependence of Ha produces an extra term, 
fi'^^dpHa). This term may be the origin of the discrepancy. If it were to 
contain the negative of the momentum entanglement entropy, it would can- 
cel the contribution of the latter to the von Neumann entropy term, leaving 
us with no momentum entanglement in (H9|) . We leave a more complete 
understanding of this point for future work. 

5.5 Non-locality of the effective action 

One further conceptual issue arising for interacting fields is that the Wilso- 
nian effective action for the lower ip modes, Sk-, is in general nonlocal. This 
raises further questions for the interpretation of the IR contribution (119|) . 
The Hamiltonian formalism corresponding to a nonlocal action is at least 
nonstandard, so the canonical thermal ensemble is nonstandard, and there- 
fore the relation between the path integral with a nonlocal action and the 
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thermal partition function is unclear. (The nonlocality problem is less severe 
at the level of the gravitational effective action F^, where a curvature expan- 
sion has the first nonlocalities appearing at order [19], beyond the range 
of our approximation.) We leave clarification of this issue also to future work. 

Finally, it is worth mentioning that in the setting of the analysis of 
Ref. [22], the entropy contributions were studied at scales above or below 
mass scales in the theory. Presumably in that setting the presence of the 
mass threshold suppressed any non-locality in the Wilsonian effective action, 
in addition to any /3 dependence of the Hamiltonian analogous to Ha ( 150|) . 

6 Summary and Discussion 

The aim of this paper has been to investigate black hole entropy within the 
framework provided by the renormalization group, in order to probe the role 
of entanglement entropy in a setting where its contribution is inherently fi- 
nite. The idea was to avoid the regulator dependence that arises for an 
otherwise divergent quantity. In Section 2, we first reviewed how the entropy 
of the canonical ensemble containing a spherical black hole is computed (ne- 
glecting metric fluctuations) from the on-shell evaluation of the full effective 
action for the gravitational field, T[g]. We then reviewed how the contribu- 
tion from a minimally coupled matter field is formally equal to its (divergent) 
entanglement entropy. 

In Section 3, we introduced an RG cutoff scale k, and defined a flowing ef- 
fective action F^ [g] for the metric, which excludes the effects of IR excitations 
that are below the scale k. The Bekenstein-Hawking entropy computed from 
this effective action, in terms of the running coupling Gk, is complemented 
by a contribution from the remaining, unintegrated modes to give the total 
entropy, according to fl27|) in general, and f l36|) for the massless scalar. In 
Section 4, we developed a similar decomposition of the entropy in the case of 
interacting quantum fields. The upshot is equation fHHj) . which differs from 
f l27|) in that the quantity encapsulates the contribution of the lower modes 
not by an explicit one-loop determinant f l2B]) . but implicitly through a path 
integral involving the Wilsonian effective action for the low-energy modes, 
Skldif] (147|) . The combination Sk = F^ + Sk evolves with k according to 
( 145|) . which is a curved space version of the Polchinski equation. The results 
of Sections 3 and 4 can in principle be extended to include gravitational fluc- 
tuations, using a background field quantization of gravity along the lines of 
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Ref. [23]. 

Section 5 was devoted to analyzing whether the contribution from the 
modes below k can be identified with their entanglement entropy as a sub- 
system. We identified a number of problems for this interpretation. Some 
just concern the precise definition of the entanglement and are thus perhaps 
not very significant, while others may pose a serious challenge to the very 
notion of an RG scale dependent entanglement entropy. Five issues were dis- 
cussed in Section 5. The first three are relevant for both free and interacting 
fields, while the last two arise only in the presence of interactions: 

(i) The well-known presence of a contact term in the entropy for non- 
minimally coupled fields means that the contribution of the lower modes 
cannot refiect only their entanglement (unless the proposal mentioned 
in footnote |9] is correct). However, the contact term is at least isolated 
from the rest of the contribution. It can be thought of as an quantum 
correction to the Noether potential at scale k, and thus as part of the 
"gravitational" entropy at that scale. 

(ii) Our use of Euclidean momentum cutoff means that the RG scale has 
no direct real space interpretation, although this may not be a serious 
impediment since a precise Lorentzian correspondence could perhaps 
be established, or the RG scale could be implemented in a different 
fashion. 

(iii) The fuzziness of the horizon concept, and therefore of the horizon en- 
tanglement entropy, in the presence of a momentum cutoff is a more 
basic issue. However this could be looked at as a necessary ambigu- 
ity in the notion of scale dependent entanglement entropy, and not a 
fundamental problem with that notion per se. 

(iv) Interactions of the matter field produce an entanglement between sub 
and super k modes that is not included in the contribution of the sub 
k modes to the total entropy of the thermal ensemble. Perhaps the 
momentum entanglement can be isolated by its coupling constant de- 
pendence. Also, if the coupling constant A is much smaller than unity, 
the momentum entanglement could perhaps be suppressed by focus- 
ing on a sufficiently small neighborhood of the horizon while remaining 
compatible with the momentum cutoff. 
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(v) After integrating out the super k modes in an interacting theory the ef- 
fective action must be non-local, since the dynamics of the sub k modes 
is not truly autonomous. This non-locality might invalidate the link 
between the sub k partition function on a cone and the entanglement 
entropy, since then no standard Hamiltonian for the system exists. 

Our conclusion is that when the gravitational black hole entropy is derived 
from the Noether charge for an effective action at scale k, the finite remaining 
contribution to the total entropy from the IR quantum modes below this scale 
has no straightforward interpretation as entanglement entropy. However, for 
free fields this interpretation may be admissible provided that difficulties 
(ii-iii) can be suitably finessed. For interacting fields the points raised in (iv- 
v) raise a larger challenge to the entanglement interpretation. In any case, 
blithe claims involving that interpretation should be avoided. 
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A Polchinski equation on curved back- 
grounds 

In this appendix we derive the RG fiow equation (l45l) that Sk[g,(t>] satisfies. 
We consider the definition (138|) evaluated for two close RG scales k and 
k + Ak, subtract the equations and expand for small Ak. obtaining: 

AA; 

Sk+Ak[9, 0] - Sk[g, 0] = — Tr ■ 7^fc ■ v?)</,,fc • (51) 

We use a compact notation, with the overdot being a /c-detivative, writing 
F ■ D = j dx F{x) D{x, y) for a function F and an operator D, and 

{F)6 k = ■' \ 7 ^ • (52) 
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The right hand side of 0511) can be related to the functional derivatives of 
Sk[g,(t>\ with respect to 0, using the following relation between two-point 
operators computed by differentiation of f p8|l : 



6Sk SSk 



Dk- {if - Dk Dk ■ (p) 



5<P6<P 50 6<P "-^"^ ^^^^ 

where stands for the two-point operator (— -|- TZk)- This leads in the 
hmit AA; — )■ to the flow equation: 



5<p 



Tr 



+ Tr 



(54) 



Using ( l43|l . this can be rewritten in a more compact way in terms of the 
low-momentum propagator P^: 



Sk[g,' 



5<p ' ^' 5(P 



Tr 



Pk- 



+ Tr 



Pk(-V 



(55) 

This is the Polchinski equation in a curved-background setting. Standard 
presentations of the Wilsonian renormalization group [lH HI] are restricted 
to flat space and include only the first two terms, omitting the third one 
since it affects only the gravitational effective action and is thus irrelevant in 
flat space. (The same thing happens for the normalization factor Nk). The 
expression (15^ in terms of the cutoff function TZk highlights the similarity 
to our framework for free fields; note in particular that if Sk[g,(p] does not 
depend on 0, as happens when the bare action is a free massless field, the 
first two terms of fl54p vanish and we recover (1251) . (The form of the equation 
is the same whether the overdot stands for or for kdk)- 
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